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ABSTRACT 

This  paper  treats  the  buckling  and  postbuckling  behavior  of  circular 
arches  and  rings  under  constant  directional  pressure.  New  exact  and 
approximate  solutions  are  given  for  the  linearized  eigenvalue  problem. 

It  is  clearly  demonstrated  that  the  assumption  of  inextensibility  is 
quite  reasonable  for  the  asymmetric  buckling  of  steeper  arches  and  of 
rings.  Asymptotic  analyses  of  early  postbuckling  behavior  are  given, 
based  on  the  theory  of  Koiter  and  the  formalism  of  Budiansky  and 
Hutchinson.  The  postbuckling  behavior  is  shown  to  be  stable,  and  un¬ 
affected  by  the  assumption  of  :.nextensibility. 
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perturbation  solution  parameters  defined  by  Eq*. (40)  and 
(37) 
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displacement  coefficient  («*  ) 

extensional  strain 

Young's  modulus 
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thickness  parameter 

moment  of  inertia 
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potential  energy  functionals 

load  per  unit  length  of  arch 

extensibility  parameter  defined  by  Eq.  (16) 

arch  radius 

nonlinear  extensional  strain 
tangential  displacement,  dimensionless 
tangential  displacement,  dimensionless 
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half  angle  of  arch  vertex 
first  order  variation 
porturbntion  parameter 
arch  coordinate 
arch  loading, dimensionless 

rotation  of  normal  to  arch  centerlines,  dimensionless 


''"in i  til  m vdL. 
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1 .  INTRODUCTION 


This  paper  treats  the  buckling  and  the  postbuckling  behavior  of 

pinned  arches,  under  constant  directional  ("djad")  pressure,  where  the 

arch  radius  of  curvature  is  a  constant.  A  linearized  treatment  of  this 

* 

problem  has  been  given  by  Chwalla  and  Kollbrunner  [l]  ,  who  assumed  that 
the  arch  was  inextensible  during  buckling.  Kammel  [2]  has  shown  that 
extensibility  did  not  have  a  great  effect  on  the  buckling  loads,  a  result 
confirmed  by  the  present  study.  Also,  Levy  [3]  and  Carrier  [4]  have 
given  "closed  form"  solutions  for  the  postbuckling  of  inextensional 
circular  rings  in  terras  of  (implicit)  elliptical  integrals.  It  is  also 
worth  noting  the  recent  contribution  of  Singer  and  Babcock — [-£r}'7  wKo^ 
discuss  the  relationship  of  semi-circular  arches  under  dead  pressure  tc 
that  of  completa  rings,  pointing  out  the  need  to  account  for  rigid  body 
displacements . 


The  present  work  treats  the  linearized  buckling  problem  both  exactly, 
within  the  framework  of  ring  theory,  accounting  for  extensibility,  and 
anoroxiiMt®"<v.  extending  levy's  approximation  for  lines  to  arches  of  various 
vertex  angles.  A  qualitative  measure  of  the  midplane  extension  is  also 
provided  from  the  new  exact  solution. 

Two  asymptotic  analyses  of  the  early  postbuckling  behavior  are  presented, 
both  based  on  the  theory  advanced  by  Koiter  [6] .  The  first  makes  use  of  the 
formalism  expounded  by  Budiansky  and  Hutchinson  [7,8,9],  and  is  carried 


*  Numbers  in  bracnets  refer  to  entries  m  the  list  of  References. 
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out  to  include  terms  in  the  displacement  expansion  of  order  e  ,  where  e 
is  the  small  buckling  deflection  parameter.  It  is  shown  that  the  post- 
buckling  behavior  of  the  arch  following  asymmetric  buckling  is  stable, 
and  it  is  not  appreciately  effected,  qualitatively  or  quantitatively,  by 
the  inextensibility  assumption.  Further,  it  is  shown  that  the  inclusion 
of  the  second  order  terms  in  the  displacements  has  no  appreciable  effect 
on  the  quantitative  outcome. 

Finally  an  anlysis  based  directly  on  Koiter's  approach  is  presented, 
also  accurate  m  the  displacement  expansion  to  order  £  .  JEhe-xesults  are. 
identical,  for  .Ihi1,  jmmIiImji;  1 1  ii iT^T~TiVrr7huMl  using  the  Buaiansky- 

Hutchinson  formalism,  and  serve  to  confirm  the  idea  that  (at  least  for 
"simple"  problems)  accurate  quantitative  information  can  be  obtained  using 
only  the  buckling  mode  shapes  in  the  ]x>stbuckiing  analysis.  Koiter  has 
previously  considered  the  elastica  [6,  and  a  two-bar  frame  [10]  using  only 
the  buckling  mode,  while  Haftka  and  Nachbar  [ll]  have  compared  results  of 
a  "one-term"  solution  with  an  exact  solution,  and  v~.he  comparison  was  quite 


good. 


”  5  " 
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a2 

— r-  +  AH 
d<T 
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,4  2 

h.  -  K  — r  +  AH  ~  +  i 

d$  d<r 


It  la  immediately  of  interest  to  nov*s  that  equations  (7)  car  be  uncoupled 
to  yield  s  single  equation  in  v,  ,  or  the  identical  equation  in  alone. 
Furthermore ,  that-sxngle  equation  is  identically  that  given  by  Chwalia  and 
Sollbrunnor[ll ,  which  they  derived  by  assuming  inextensxbility  in  the  buckling 
state,  i,e. ,  e^  *  0.  That  equation  is,  obtained  here  with  no  assumption 
regarding  the  extensibility? 


<LaV  V 


Vi 

V1  T 


(2  +  A)v 


XV 


tt  +  2.\)v^' 


+  Av, 


(9) 


The  prime  superscripts  in  equation  (9)  correspond  to  differentiation  with 
respect  to  the  arch  coordinate.  The  question  of  inextensibility  will  be 
raised  again  in  the  sequel. 

For  completeness  the  natural  {in  terms  of  stress  resultants)  and  geo¬ 
metric  boundary  conditions  for  the  linearized  prcolera  are  also  noted; 


Either 

N!  -  M1  "  0 
dM 

or 

V1 

is  prescribed 

Either 

5T "  **i  =  0 

or 

W1 

is  prescribed 

dw.. 

Either 

Mi =  0 

or 

L 

is  prescribed 

The  boundary  conditions  are  applied  at  <J>  ~  i  a,  where  2a  is  the 
arch  vertex  angle  {Figure  1).  Results  for  the  complete  ring  are  ol  tam¬ 
ed  by  setting  a  -  ~  [5], 


3.  LIMEARI2ED  (BUCKLING)  SOLUTIONS 


A  solution  to  the  eigenvalue  problem  represented  by  the  differential 
equation  (-3)  has  been  given  by  Chwaila  and  Kollbrunner  [1]  for  a  simply 
supported  arch,  using  the  boundary  conditions  {reflecting  the  inextensibility 
condition) 

«  III 

V1  =  V1  =  VJ  =  0  at  4  =  ±  a  (11) 

Some  of  the  resulting  eigenvalues,  corresponding  to  a  solution  that  is 
asymmetric  with  respect  to  $  =  0,  are  shown  in  Table  1. 
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(14) 


a 

Some  numerical  values  of  A  are  given  in  Table  1.  The  agreement  is 

rr 

obviously  quite  good.  Also,  At  is  noted  here  that  for  the  case  v  =  i/2, 
i.e.,  a  complete  ring,  the  mode  shape  represented  by  Eas.  (12)  has  been 
shown  to  be  essentially  identical  to  that  obtained  by  Chwalla  and  Kollbrunner 
in  their  analysis  [5] . 


Finally,  in  an  attempt  to  assess  the  effect  of  the  inextensibility 

•• 

assumption  on  the  buckling  load  (see  the  work  of  Kammel  [2])  and  to  lay  the 
groundwork  for  the  study  of  the  effects  of  inextensibility  on  the  post- 
buckling  behavior,  an  e«cct  solution  to  the  differential  equations  (7)  ana 

the  boundary  conditions  (10)  is  now  given.  The  asymmetric  solution  to 

* 

equations  (7)  is 


W1  =  A1  s*n  I7«  +  sin$  +  cos$ 


(15) 


t, 

V1  = 


H1 

-  ■=.  COS  -/a*  $  +  (QA_  -  A _)cos$  +  A  <t>  sin <f> 

'Tv  J 


where 


q  =  LUiiL. Z...U.  (16) 

C  i  -  H(A  -  1) 

and  the  A^  are  constants,  determined  to  within  an  arbitrary  factor  by 
satisfaction  of  boundary  conditions. 


*  This  solution  is  not  valid  for  A  =  1.  However  for  O  <_  $  <  ~/2,  thi-  is 

of  no  concern  here,  as  A  >  3. 


8 


Using  the  results  (IS) ,  the  exact  in-plane  strain  and  rotation  formulae 
(for  buckling)  can  be  written  as: 


e1  =  (1  ~  Q)  A3  sin$ 


xf  =  ('ft1-  —^AjCOsJT^  t  (1  +  Q)A^cos$ 


(17a) 


(17b) 


The  eigenvalues  (buckling  loads)  are  obtained  after  satisfying  the  following 
boundary  conditions: 


EE  E  A1 

wi  “  vi  * 


o  at  $  «  *  a. 


The  resulting  transcendental  equation  from  which  the  buckling  load  is  calculated 


is  found  to  be: 


'  *  (1  +  2Q  -  AQJsinacosa  +  a(l  - 


:(1  -  A) 


For  the  sake  of  comparison,  the  eigenvalue  equation  obtained  by  Chwalla  and 
Kollbrunner  is  also  given. 

4a*  tan  a  = - - -  (20) 

(3  -  A)sinacosa  +  <x(l  -  A) 

It  is  seen  that  Ects.  (19)  and  (20)  are  identical  when  Q  =  1,  which  corresponds 
to  inextensional  deformation  (Eq .  (17a)).  Numerical  values  for  the  buckling 
pressures  given  by  Eq.  (19)  are  presented  in  Table  1,  and  again  the  agreement 
is  excellent. 


It  is  worth  noting  at  this  juncture  a  few  features  of  the  extensional 
strain  e^  resulting  from  the  exact  solution.  First,  considering  the  co¬ 
efficient  on  the  right  hand  side  of  Eq.  (17a),  and  the  definition  of  Q  (Eq. 
(Ib^),  it  can  be  seen  that 


r 


-  9  - 


i  -  Q  = 


2H(A  ~  1? 

1  -  H(.\  -  1) 


H  «  1. 


(21) 


proportional  to  the  square  of  the  thickness-fco-ra&ius  of  the  arch,  and 
is  thus  quite  snail. 

Secondly,  it  can  be  observed  fron  Ea.  (17a)  that  the  exfcensional 
strain  is  an  odd  function  of  the  arch  coordinate,  i.e., 

e*  {-  $)  =  -  e®{$)  (22) 


Also,  it  turns  out  that  there  exists  another  uncoupling  of  the  exact 
differential  equations  (7)  that  yields  the  result 

,2  E 

d  E  .... 

— 2 —  +  e^  =  0  (2.»; 

Thus,  in  addition  to  the  odd  function  represented  by  the  sinusoid  al¬ 
ready  obtained,  there  is  another  solution  to  the  equations  for  which 
£ 

w^(d)  and  e^{$)  are  even  functions  of  $  ,  tl  &  letter  being  a  cosinusoid. 


This  solution  is  not  considered  at  all  in  the  present  work  as  attention  is 
restricted  oriy  to  arches  steep  enough  for  as/metric  (bifurcation)  buckling 
to  take  place. 
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4.  POSTBUCKLIliG:  TIE  BUDIANSKT-HUTCHINSON  APPROACH 


In  this  section  the  early  postbuckling  behavior  of  the  ar^h  will  be 
obtained,  following  the  asymptotic  analysis  outlined  in  the  very  lucid 
exposition  of  Bv 'iansky  [9j.  As  with  the  buckling  analysis,  only  asym- 
metrij  buckling  of  steeper  arches  is  considered  here.  For  very  shallow 
arches  symmetric  buckling  can  take  place  at  lower  pressures  [15}.  how¬ 
ever  this  snap  type  of  buckling  cannot  be  analysed  with  the  Koiter  theory. 
The  starting  point  is  the  virtual  work  (or  variational)  statement  of  the 
full  nonlinear  problem: 
a 

/  [N  $e  +  NX  6x  *  M  6k  -  Xdw]d$  =  O  (24) 

-a 

In  addition  to  the  nomenclature  defined  in  the  previous  section,  the 
curvature  K  has  been  introduced  in  2a.  (24) .  Then 


M  =  K  =>  - 


d£ 

d  $ 


(.15) 


The  variational  equation  (24)  must  hold  fcr  any  kinematically  consistent 
set  of  displacement  variations  6v  and  6w.  with  corresponding  strain 
field  variations  £e,  dx  and  5k. 


The  following  expansions  are  now  introduced: 


v  »  Xv  tv,  +  e  v. 

o  1  2 

2 

w  *  Xw  +  ew,  +  e  w„ 

o  1  2 

2 

e  =  Xe  +  ee,  +  e  e0 

o  1  2 

X  =  XxQ  +  exx  +  c2x2 


K  = 


CK,  t  £  K_ 
1  2 


^  •  •  •  • 

+  .... 

*  .  .  * 

+ 

f 


.... 
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N  =  -  K  +  eN,  +  e2N„  +  . . . . 

O  l  2 

H  =  eMx  +  e2 M2  (26) 


Here  c  is  a  small  parameter  such  that  when  e  -*•  0,  then  X  •+  X  , 
the  biturcation  buckling  pressure.  The  curvature  and  moment  expansions 
begin  with  terms  of  order  e  to  denote  that  the  pre-buckling  state 
(zeroth  order  terms  in  Eqs.  (26))  is  a  membrane  staie.  In  fact,  as  it 
is  easily  demonstrated  that  a  nonlinear  membrane  state  cannot  exist  for 
this  problem,  the  quartity  XQ  shall  be  henceforth  taken  as  zero,  and 
thus  only  a  linear  membrane  pra-buckling  state  is  considered. 

For  the  "unvariedr  quantities  in  Eq.  (24)  the  expansions  (26)  are 
substituted.  After  rearrangement  by  increasing  powers  of  e  ,  the  following 
new  variational  equation  results: 


/  [-  tl  5e  -  A6w  ]d$ 

-a  ° 


a 

+  t  f  [N,  6e  +  M,<5k  -  N  X,  $x]d$ 

_  1  i  Ox 

— Q 


+  /  rtf 

-a  *  2  ~ 


M  Av  x  £vr  v  B  vt  w  • 

•*2  ®  *,qA2> 


(27) 


3  a 

+  e  /  [N3oe  +  «36k  +  (N^  +  -  N0X3>  6xjd<}» 


=  0. 


The  pre-buckling  state  is  defined  by  the  vanishing  of  the  zeroth 
order  variation.  After  a  simple  calculation  that  state  is  found  to  be 
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N  a  x,  e  =  -  H,  w  =  H,  V  =  0  (25) 

o  o  o  o 

With  the  vanishing  of  the  zeroth  order  terra  in  Eq.  (27) ,  an  e  may 
be  divided  through,  and  thus  the  buckling  state  itself  is  now  found  to 
be,  recalling  that  e  -*■  0  implies  X  -*■  A  : 
a 

/  &J.  <$e  +  M.  OK  -  X  X,  6x]d$  «  O  (29) 

_  l  1  cr  l 

-a 

The  Euler-Lagrange  equations  of  this  variational  equation  are: 


dM, 

*crXl  - 

O 

n 

IS 

(30a) 

dX, 

d2M. 

cr  d$ 

*  2  *  0 
dy 

(30b) 

Expressed  in  terms  of  displacements,  Eqs.  (30)  are  identical  to  the 
linearized  equations  (7)  given  earlier. 


It  was  pointed  out  earlier  that  the  set  of  displacement  variations 

constantly  alluded  to  need  only  be  kinematically  consistent,  and  since 

the  buckling  displacements  most  surely  are  such  a  set,  they  can  be 

appropriately  inserted  ir.tc  Eq.  -29)  to  yield  the  Rayleigh  Quotient 

(Budiansky's  "energy  equation")  for  the  critical  load, 
a 


X 

cr 


/(«1e1  + 

■c _ 

a  2 
/  (XL)  d$ 


(31) 


The  same  substitution  in  the  more  complete  variational  statement  (27) 
must  be  equally  valid,  and  noting  the  quotient  (31) ,  it  is  found  that 


l' Ue?*v3'"-»T' 
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(X  -  X)  /  xt  <1$ 
v  cr  '1 

~a 


+  f.  /  [ti^  +■  M2Kx  +  Kx,l  -  XxiX?)d« 


2  K  7 

+  c  f [«3ex  +  M3Kj  +  N.x1X2  +  H2x1  '  A>‘Ix^d*  =  0 

-  a 


At  this  point  the  relationship  between  (X  -  /  appears  to  be  taking 

shape.  However,  it  appears  that  equations  of  the  third  order  will  have  to 

2 

be  solved  in  order  to  obtain  the  r oefficient  of  e  above.  In  order  to 

* 

avoid  this  difficulty,  an  orthogonality  condition  is  introduced  as  follows 


/  [ll e  +  M.K  ]d«  «  0;  n=2,3,4. 

_  l  in  jl  m 
-a 


It  follows  from  Eqs,  (33)  and  (29) ,  since  the  m  order  mods  shapes  must 
be  kinematically  admissible,  that 


r*iV  ■  °!  n  *  2-3‘‘ 


in  order  to  apply  the  orthogonality  condition  i33'  and  xts  coiisequiutc 
(34)  successfully,  the  following  consequences  of  constitutive  linearity 
and  geometric  nonlinearity  must  be  noted: 


*  The  precise  meaning  of  this  orthogonality  condition  has  not  been 
established  and  is  the  subject  of  some  discussion?  see  Reissner 

[14]. 


9*‘  ~g*_,  _  ^  _  _ .  „ 
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N2el  = 

N2S1  * 

N1S2 

+  2  Xl); 

N3S1  “ 

N3Si  = 

"A 

'  Ve3 

+  XjXj) i 

(35) 

M2Ki  “  M1K2?  M3K1  =  MlK3f 

In  Eqs.  (35),  S^  has  been  used  to  denote  the  total  nonlinear  strain  of 
order  e1  ,  which  is  made  up  in  part  of  lower  order  terms  due  to  the 
geometric  nonlinearity.  Then  notings  Eqs.  (35),  (33)  and  (34),  the 
variational  equation  (32)  car.  be  put  into  the  very  simple  form: 


Since 


X  = 


\cr  :•  «•»+!>■€'■+... 


(36) 


where 


a,  „  ! 
3 


/  N  xf  d* 

-tt  -1 

a  2 

/  X7  d* 

-a 


(37a) 


a 

/ 

-a 


/  [->niXix2  +  Vi  ]d* 


a  2 

/  xf  d$ 

-a 


(37b) 


„  1  _  1  .  „  1  _  _  1 ,  .  1  2 

H1  "  H  ~1  H  =1  N2  ~  H  S2  “  H(e2  2  Xl} 


(38) 


the  result  (36)  can  be  written  in  its  final  form: 


X_ 

X 


«  1  +  ae  + 


be2  + 


cr 


where 


a  =* 


a' 


cr 


2BX 


cr 


/elX,  d$ 

tr 

-a 


(39) 


(40a) 


(40b) 


b  = 


/72elXlX2  +  e 2 ^ 


— 

~  f  i  J®-" 


cr 


HA 


cr 


d* 


-i 


Recall  that  according  to  the  theory,  it  a  *  G,  then  the  structure  will 
be  in  an  unstable  state  in  the  poi.tbucM.ng  range,  and  it  will  also  be  im¬ 
perfection  sensitive.  If  a.  -  0,  \  f  0,  then  the  postbucklmg  stability 

and  the  imperfection  sensitivity  vary  as  the  arithmetic  of  b. 


In  the  present  result  it  appco~s  at  first  glance  that  the  existence  of 
a  stable  or  unstable  equilibrium  state  in  the  pcstbucklmg  range  depends  on 
whether  or  not  inextensicility  is  assumed.  For  from  equation  (40a ) ,  if 
e^f^)  is  non-zero,  then  the  coefficient  a  might  be.  However,  recalling 
Eq.  (22;,  it  is  clear  that  the  coefficient  m  question  will  alwavs  vanish 
for  asymmetric  buckling !  Thus  the  exclusion  (or  inclusion)  of  the  arch 
extensibility  can  at  best  only  influence  the  quantitative  nature  of  the  post- 
buckling  equilibrium  state,  and  not  its  qualitative  nature 

It  has  already  been  shown  that  for  asymmetric  buckl.ng  the  buckling 
extensibility  is  quite  small  in  magnitude,  with  no  appreciable  effect 

on  the  buckling  load.  With  this  in  mind,  as  well  as  the  discussion  of  the 
previous  paragraph,  it  seems  quite  reasonable  to  stipulate  that 


e^  -  0  and  e^  =  O 


In  this  instance  the  coefficient  b  reduces  to 

f  X?  d$ 

b«  1 


2H>. 


cr 


f  xj  dj 

-a 


(41) 


(42) 
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Using  the  excellent  approximate  solution  (12) ,  with  its  critical  load 
defined  b_  Eq.  (14) ,  it  is  then  a  straightforward  matter  to  calculate 


the  postbuckling  relationship: 


A  =  A  +  (12)  (~) 
cr  z 

it 


[1  +  3  (4-d + §<4  -  D2'] 

_ o _ a _ 

[2  +(4  -  l)2] 

a 


* 

(£A_)2 
4  h  ’ 


(In  writing  this  result  cognizance  has  been  taken  of  the  fact  that  the  co¬ 
efficient  A  in  Eq.  (12)  was  nondimens ionalized  with  respect  to  the  radius 
* 

A 

R,  i.e.  a=  —  ,  and  of  the  definition  (3)). 

R 


It  is  clear  from  Eq.  (43>  that  A  >_  A  #  and  thus  the  postbuckling 
behavior  is  always  stable,  and  that  the  arch  (or  ring  if  a  =  will  not 
be  sensitive  to  imperfections  :in  its- behavior. 
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5.  POSTBUCKLING ;  KOITER'S  APPROACH 

In  order  to  examine  the  stability  of  the  bifurcation  point  and  of 
the  postbuckling  behavior  using  Koiter's  approach  directly,  the  total 
potential  energy  is  first  divided  into  the  prebuckling  contributions, 
and  into  the  terms  due  to  the  buckling  and  prebuckling  behavior.  If 
the  prebuckling  behavior  is  again  linearized,  the  "expansions"  (44)  are 
now  introduced,  i.e.,  let 

e  =  Ae  +  e 

o 

w  =  Aw  +  w  (44) 

o 

X  »  X 

Substituting  these  expansions  into  the  energy  functional  (5)  yields 

f  »  +  V)  -  2A2wo)d4» 

-a 

1  ° 

x  i.  r  /oi -  m 

t  ??  J  Q  *  4An«v/uy 

H  O 

-a 

,  1  a.~  r  2  .  2  2  1  4,  . . 

+  /  [e  +  AeQX  +  ex  +  j  X 

-a 

+  /  (|£)2  d<J>  (45) 

-a  9 

In  view  of  the  prebuckling  results  obtained  earlier,  the  total  potential 
energy  due  to  deformation  beyond  the  prebuckling  state  can  be  written  as 
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||  (U  +  V)  +  2aHA2  =  P*  +  P*  +  P*  +  P* 
El  12  3  4 


(46) 


where 


Px  »  -  2A  /  (e  +  w)d$ 

-a 

—A  1  r  2,,  -  “  2, , 

P2  =  K  f  *  ~  *  f  X  &<t> 

-a  -a 


/&>2 


-a 


1  7  2., 

-  /  ex  d$ 

-a 


3 

A  _  1  7  4 

P4  "  4H  X 
-a 


Now,  by  the  definition  of  e($)  it  is  clear  that 
c 

c  -  2x  /  (^  -  w  +  w)d$  =  -  2>jv(a)  -  v(-a)j 


(47a) 

(47b) 

(47c) 

(47d) 


(48) 


A  _ 


Hence,  for  a  supported  arch,  P^  =  0.  Thus  the  potential  energy  change 
during  and  after  buckling  is  simply 


(49) 


PX  «  Pi;  +  P*  pA 

2  3  4 

where  Px,  px  and  ?x  are  given  by  equations  (47b,  c,d) . 

Recall  that  the  behavior  of  interest  takes  place  at  loads  very  close  to 


A  _  A 


cr' 


t*ho  KnrVl  inn  lnarl 


-»  * - 


Htution  t**c  riisr^y 


functionals  (47)  can  be  expanded  in  Taylor  series  about  the  point  X  t-.  X 
X  X 

As  P^,  P^  do  not  (for  this  problem)  explicitly  depend  on  x,  they  appear 


cr 


only  as  constants,  while 


P2 ' 4=!  +  <»  -  Kr) 

'  cr 


dP. 


dA 


A-A 


cr 


'V“-  XcrIP2  +  - 


S50) 
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where 


a  a  a 

P,  =  ~  /  e2d4>-  X  /  X2d<}>  +  /  M)2*$ 

2  H  -a  cr  -a  -a  d* 


(51a) 


Pi  =  -  /  X  d* 

-a 

Then  the  total  energy  change  is 


(51b) 


P  -  P2  +  (X  -  Acr)  P'  *  P3  +  P4  (52) 

where  P_  =  P^  and  P.  =  P*  . 

3  3  4  4 

Now  the  following  displacement  and  strain  expansions  are  introduced: 
o 

w  =  ew^  +  e*w2  +  • • • 

2 

v  =  ev  +  e  v  +  ... 

1  2  (53) 

e  =  ee^  +  e2e2  +  ... 

X  =  +  e2X2  +  .  •  . 

where  £  is  the  same  small  parameter  used  previously.  Introduction  of  the 

4 

expansions  (53)  into  the  energy  functionals  yields,  to  order  e  : 


a  a  a  dx 

p2  *  (i  1  eld$  '  Xcr  1  xld*  +  !  (ar)2  ddi 

-a  -a  -a 


+  2G  {H  X  eie2d*  "  Xcr  f  X1  X2d*  +  1 


a  dXjL  dx2 


d4>  d  $ 


(54a) 


4  1  0  ?  °  2 

+  e  %  f  'e2  +  2eie3,d(J’  **  Xv  ^  <X2  +  2X^)5^ 
-a  -a 

a  dx0  _  dXi  dx, 

+  +2dT5r’d*> 
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«  ,  a  ,  4a2 

P2  =  "  e  1  Hd*’  2e  f  XiX2d'*»  -  e  /[x2  +  2x1X3J<S4' 

-a  -a  -ci 


(54b) 


P„  = 


e~  H  eiXia^  +  e  h  1  te2Xl  +  2eiXiX23dl^ 
-a  -a 


(54c) 


P„  = 


41  ;  4_. 

E  4H  '  *ia* 
-a 


(b4d) 


Then  the  same  orthogonality  conditions  applied  in  the  Budiansky- 
Hutchinson  formalism  are  applied  again  (written  entirely  in  terras  of 
kinematic  quantities)  to  tne  functionals  (54) .  These  functionals  are  then 
greatly  simplified  and  take  the  form: 


P 


2 


P' 

2 


=  e2(  i  /  e2d*  -  X  /  Sc2d*  +  f  <^)2d*> 
.  ,  o.  a  „  a  dx_  ,, 

+  e4{«  1  cld*  ~  V  f  x?d*  +  f 

H  -a  2  cr  -a  2  -a  d* 

2  r  2aa  4  °l  2,, 

=  -  e  /  X^  -  e  /  x2^<J 

-a  -a 


(55a) 


(55b) 


p3  -  e3  £  /  exXi  d« 


P.  = 


+  e4  £  'le2Xx  +  2elXlX2]d^ 

**0t 


4  1  4 

c  h  1  h  d* 

-a 


155c) 


C55d) 


2 

Note  that  the  terra  of  order  e  in  P^  is  simply  the  Rayleigh 

quotient,  and  thus  it  vanishes  identically.  Further,  the  terra  of  order 
4  ’ 

e  in  P2  will  be  multiplied  by  (A  - 


according  to  equation  (52) . 


Since  it  i»  expedited  that  for  this  problem  that  (X  ~  X  )  will  be  of 
2  4  * 

order  £  ,  the  term  of  order  e  in  ?2  may  be  deleted  as  being  of 
higher  order.  Thus  the  fina1  form  of  the  energy  change  is: 

a  a 

P  «  -  e2{X  -  Xcr)  /  +  c3  |  /  elXJ  d$ 

-a  -a 

a 

+  e4  {i-  /  [e2x^  +  2e1X1X2  +  |  X^]  d$  (56) 

—a  *" 

♦  i  /  -  A  /  X2d«  +  “  } 

—a  -a  ~o 

The  resemblance  between  the  present  approach  and  the  Budiansky-Kutchin3on 
approach  is  now  more  visible.  However,  it  appears  that  there  will  be  some 
difference  due  to  the  term 

a  o  a  dX 

P2  [w2,v2]  -  |  /  e2d*  -  Xcr  /  +  /  (^>2d^  (57) 

—a  -a  -a 

4 

in  the  expression  of  order  e  in  Eq.  (56) .  Terms  o*  this  type  do  appear 
consistently  in  Koiter's  original  analysis  (see,  for  example,  Section  3 *»  of 
Ref.  [6].),  although  they  do  not  appear  in  analyses  using  the  Eudiansky- 
Hutchinson  formalism  (see,  for  example.  Ref. [9].).  For  the  present  work, 
under  the  well  established  assumption  of  inextensibility,  the  terms 
P2^w2'v2^  ke  shown  to  yield  no  contribution  to  the  result,  and  the 

results  will  be  identical  to  those  already  obtained. 


The  cecond  order  displacements,  v2  and  w2,  are  obtained  as  the 


solution  of  the  variational  aquation 


6(i  '4a-  +  fi ^r>2« 

-a  -a  -a 


+  ~  /  le2X^  +  2e1y1x2  +  A  j^d*}  *  O 


In  carrying  out  this  variation  all  terns  of  first  order  are  of  course 
considered  as  known  or  given  quantities.  If  inextensibility  of  the  arch 
is  assumed,  ei  "  e2  ”  then  the  variational  equation  is  just 


“  dX?  2 

«{/ 


*  /  X?d<-}  =  0 

C*.  _  c 


1'hus  the  variational  equation  leads  to  set  of  homogeneous  equations  that 

* 

can  have  only  a  trivial  solution  .  hence ,  for  the  inextcnsible  arch,  it 
is  passible  to  take  as  solutions  for  the  second  order  displacements  the 
trivial  ones,  i.e,, 

w2  ~  v2  =  °  (60) 

Hence  the  energy  functional  reduces  to  (recalling  -  O  again) : 

P?<  »  -  c2(X  -  X  >  /  >;2  d*  /  xjett  (61) 

-a  -a 

* 

Using  the  approximate  solution  (12) ,  with  A  =  A  /R,  it  is  a  rt.aiyht- 

forvard  matter  to  vcMpute  the  energy  change  t*  as 

*32  *52  2 

P  *  (— )  -12  ♦  -  1)  ]Ucr  -  H  *  <—>  %<-  ~1)  et-y-mi] 

v  a  £7  il  a  a 


*  A  solution  linearly  dependent  on  the  buckling  solution  is  ruled  out  by 
the  orthogonality  condition. 
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» 

At  the  critical  lead,  l  =  >  ,  it  is  sesn  chit  E  0  and  thus  »>« 

cr 

bifurcation  (critical)  state  is  itself  stable.  The  initial  postbuckling 

"k 

load-deflection  relationship  is  determined  by  the  condit*  >n 


\ 

dP  ’ 


d(eA*) 


-  t 


(63) 


from  which  it  follows  tii3t 


2  n  +  - 1>  -  f<4“ i)2]  * 

X  =  Xcr  +  (12,  % - ^ - - 1 -  (^2 

IT  V  2 

2  +  (~ - ir 

a 


(64) 


This  is  exactly  the  result  (43)  obtained  previously. 


*  This  is  an  approximation  to 


.  k 
£P 


=  O 
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6.  DISCUSSION  AND  CONCLUSIONS 


The  present  paper  has  considered  the  buckling  and  postbuckling 
behavior  of  circular  arches  (and  rings)  under  constant  directional 
loading.  For  arches  steep  enough  so  that  asymmetric  bifurcation 
buckling  is  critical,  it  has  been  demonstrated  that  the  arch  inextensi¬ 
bility  is  a  negligible  quantity  in  magnitude,  in  its  effect  on  the  buckling 
load,  and  in  its  effect  on  initial  postbuckling  behavior.  In  addition  an 
approximate  solution  for  the  linearized  eigenvalue  problem  was  presented 
that  duplicated  almost  precisely  the  buckling  pressures  obtained  from  the 
exact  solution. 

The  numerical  results  for  the  eigenvalue  problem  as  solved  by  chwalla 
and  Kollbrunner  and  as  given  by  the  present  approximate  and  exact  solutions 
have  been  given  in  Table  1.  It  is  worth  pointing  out  that  the  approximate 
(in extensible)  eigenvalues  are  always  greater  than  or  equal  to  Chwalla  and 
Kollbrunner 1 s  (inextensible)  results,  in  accord  with  Rayleigh's  Principle. 
The  results  when  the  extensibility  is  accounted  for  are  always  less  than  the 
inextensible  results.  This  is  in  direct  contrast  to  the  analysis  of  Kammel 
[2]  which  indicated  that  the  buckling  pressures  would  increase  if  extensi- 

•k 

bility  were  accounted  for  ,  As  it  has  already  been  pointed  out,  the 

*  Singer  and  Babcock  [5]  have  already  pointed  out  that  Kamnel's  result  was 

also  erroneous  in  arriving  at  the  coefficient  A  =  4  for  constant 

cr 

directional  pressure,  assuming  inextensible  deformations. 
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differential  equation  solved  by  Chwalla  and  Kollbrunner,  assuming  in¬ 
extensibility,  is  identical  to  that  solved  here,  without  that  assumption. 

The  difference  is  in  the  boundary  conditions,  where  the  inextensibility 
assumption  provides  less  freedom  {greater  stiffness) .  Hence  it  is  quite 
proper  that  the  inextensible  eigenvalues  be  upper  bounds  to  the  (freer) 
extensible  eigenvalues. 

It  is  also  worth  pointing  out  that  the  results  given  in  Table  1  for 
the  Chwalla-Kollbrunner  eigenvalues  are  not  precisely  those  given  in  the 
original  paper  [1].  Rather  they  were  calculated  (using  a  digital  computer) 
directly  from  the  transcendental  equation  (20) .  There  arc  some  dis¬ 
crepancies,  although  they  are  very  small.  For  example: 

Chwalla  and  Kollbrunner  [1]  Chwalla  and  Kollbrunner 

_ _  Eg,  (20)  -  Digital  Computer 

a  *  j  8.725  8.727 

a  =  j  3.265  3.271 

It  is  easily  verified  that  the  original  Chwalla  and  Kollbrunner  results  satisfy 
Eq.  (20)  to  within  the  precision  of  limited-accuracy  trigonometric  tables  and 
slide  rule  accuracy  that  must  have  governed  calculations  in  the  early  1930's. 

The  postbuckling  behavior  for  the  constant-directionally  arches  and  rings 
has  been  shown  to  be  quite  stable,  and  thus  imperfection  sensitivity  is  not 
an  issue  here.  Also  it  has  been  shown,  under  the  reasonable  assumption  of  in¬ 
extensible  buckling  and  postbuckling  deformation,  that  the  second  order  dis¬ 
placement  solution  does  not  change  the  results  that  would  be  obtained  if 
only  first  order  (buckling)  displacements  were  used  in  the  calculations. 
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Finally  it  has  .been  shown  that,  ror  the  special  circumstances  of 
this  problem,  the  approach  of  Koiter  as  originally  given  and  the  variation 
proposed  by  Budiansky  and  Hutchinson  lead  to  identical  results. 
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TABLE  1  ;  ASYMMETRIC  BUCKLING  PRESSURES  FOR  CIRCULAR  ARCHES;  X  «  qR3/EI. 


Solution  of  chwalla 
and  Kollbrunner 

Eq.  (20) 

Present 
Approximation 
Eq.  (14) 

Present  Exact 

Result  Eq.  (19) 

Ch/R  =  1/10) 

K 

10 

99.979598 

99.979598 

99.979541 

n 

8 

63.967766 

63.967766 

63.967676 

it 

6 

35.941318 

35.941320 

35.941151 

it 

4 

15.859006 

15.859031 

15.858590 

it 

2 

3.271245 

3.272727 

3.269233 

